Tunneling Spectroscopy of the Underdoped High-T c Superconductors 
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Considering a possibility that the pseudogap state in the underdoped high-T c superconductors 
is not due to a superconducting correlation, we study the tunneling spectroscopy below T c based 
on the phenomenological models of the pseudogap state. Specifically we consider two cases: (1) the 
pseudogap is a simple suppression of the density of states with an unknown origin; (2) the pseudogap 
state is due to an antiferromagnetic correlation. For both cases we calculate dl/dV using the BTK 
theory. The results are discussed in comparison with experiments. 
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I. INTRODUCTION 

The suppression of density of states around the Fermi 
level in the underdoped high-T c superconductors (HTSC) 
above T c is observed in various normal state experiments 
such as optical conductivity, dc-resistivity, angle resolved 
photoemission, NMR, tunneling spectroscopy, neutron 
scattering, specific heat, Raman spectroscopy etc jjj, and 
this phenomena is termed as pseudogap (PG). Currently 
a consensus of the origin and nature of this PG is still 
lacking and undoubtedly it is a key issue to be resolved 
to make any progress toward a theory of high Tc super- 
conductivity. 

About the origin of the PG state only two possibilities 
are logically allowed. The first possibility is that the PG 
is somehow related with the superconducting correlation 
and it develops into a real superconducting gap below Tc. 
And the second possibility is that the PG has nothing to 
do with a superconducting gap but with something else. 
Along the first line of thinking, preformed pair scenario 
pL pairing fluctuations scenario ||, etc. are proposed. 
For the second possibility, antiferromagnetic correlation 
|| , charge stripes H , etc are considered as the origin of 
the PG. At present experimental evidences exist both for 
superconducting origin || and for non-superconducting 
origin ^) . Therefore it is necessary to design an exper- 
iment to identify distinct features among the proposed 
scenarios. Recently we proposed an experimental test us- 
ing tunneling spectroscopy in the PG state, specifically, 
for the preformed pair scenario |lC| . Namely we claim 
that there should be an Andreev reflection signal even 
above Tc but below T* (PG cross-over temperature) if 
there exist preformed Cooper pairs without long range 
pha se coherence. Until now there are only one positive 
plf and one negative jL2j experiment reported on the ob- 
servation of an Andreev signal in the PG region of the 
underdoped HTSC compounds. 

In this paper we examine the second possibility for the 



PG state, i.e., that the PG is irrelevant from the SC 
gap. We calculate the tunneling conductance (dl/dV) at 
zero temperature when the PG coexists with a SC gap. 
Specifically we considered two cases: (1) the PG is a sim- 
ple suppression of density of states of an unknown origin 
and the SC correlation has no direct interplay with the 
PG; (2) the PG is caused by an AFM correlation (ap- 
proximated by SDW) and below Tc the AFM correlation 
and the SC correlation coexist and interplay with each 
other. For both cases, we use the BTK theory [jl3| to cal- 
culate the tunneling conductance (dl/dV). The results 
show characteristic features of the tunneling density of 
states in each case and those features can be used to sort 
out the true origin of the PG in comparison with exper- 
iments. 

The main results are: (1) for the first case, the ba- 
sic line shape of the tunneling conductance is a simple 
superposition of a standard BTK conductance and an 
assumed PG density of states. When the SC gap size is 
smaller than the PG size, the SC gap feature shows up as 
a distinguishable peak inside the PG. However if the size 
of SC gap is larger than the PG, the PG feature is over- 
whelmed by the SC gap feature; (2) for the second case, 
there is an interesting interplay between the AFM and SC 
correlations. Irrespective of the relative sizes of the PG 
and SC gap, the tunneling conductance shows only one 
gap feature at E = A to t a i (= \J A| c + Ap G ). However 
depending on the relative sizes of the PG and SC gap, 
the line shape of dl/dV looks very different. When the 
SC gap is bigger than the PG, it looks more like a con- 
ventional NIS (Normal metal-Insulator-Superconductor) 
junction. But for the other case the main feature of the 
dl/dV curve is determined by the SDW correlation and 
shows no diverging density of states approaching the gap 
energy in contrast to the NIS junction. This difference 
comes from the fact that the tunneling density of states of 
the NISDW (Normal metal-Insulator-SDW state) junc- 
tion obtained by the BTK theory (even at large barrier 
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potential (Z) limit) is not the same as the actual den- 
sity of states, which would have been obtained by the 
tunneling Hamiltonian method |Q. The difference of 
the tunneling conductance by the BTK theory and by 
the tunneling Hamiltonian method in the SDW state, 
which is in contrast to the NIS (Normal metal-Insulator- 
Superconducor) junction, shows that the description of 
tunneling process by these two methods are not the same. 
The important and interesting questions are then which 
description is more proper description for actual tunnel- 
ing process and why in the case of NIS junction those 
two methods seem to give consistent results for a Z > > 1 
limit. More details will be discussed in the later section. 



II. PG MODEL I: SIMPLE SUPPRESSION OF 
DENSITY OF STATES 

In this section we consider the case that the PG is a 
simple suppression of density of states, which has been 
developed above T c by an unknown origin. An assump- 
tion is that the SC gap develops below T c on top of the 
already developed PG density of states and there is no 
other correlation between the SC gap and the PG. 

The tunneling conductance of NIS junction is calcu- 
lated by the BTK theory as follows. 



dI(eV)/dV = 2ev F A / dE N PG (E) 



df(E- V) 



dV 



■[l + \A(E)\ 2 -\R(E)\% 



where 



\A 2 = 



\R\ 2 = 



(E+(1+2Z 2 )JE 2 -A 2 ) 2 
AS 



iZ 2 (E 2 - A 2 S )(Z 2 + 1) 
(E+(l+2Z 2 )y/E 2 ~A%) 2 
iZ 2 (Z 2 + l)(A%~E 2 ) 
^E 2 ) 



for E> A s 
for E < A s 

for E > A s 
for E < A s 



(1) 



(2) 



(3) 



are the coefficients for the Andreev and normal reflec- 
tions, respectively, and Z is the strength of the insulating 
barrier (Z = '", v ). The above equation is a standard 



formula for the tunneling conductance by the BTK the- 
ory with one modification, i.e., multiplied by the PG den- 
sity of states N PG (E). We simulate N PG (E) by Dynes' 
formula with the SC gap replaced by the PG (Ap). 



Npc(E) = 2ttN(0)R ( 



uj + iT 



A 2 



(4) 



The numerical calculation has done for a tunneling 
into (1,0,0) direction toward the HTSC and assume that 
Ap and As has the maximum value at the same direc- 
tion (1,0,0). To simulate the surface roughness we also 



present the angle averaged results by simply replacing 
As,p => Ag P cos(26>). All the presented tunneling con- 
ductance is normalized by the normal-state resistance 
R N = (1 + Z 2 )/(2N(0)e 2 v F A) and the bias voltage is 
measured in unit of Ap. Also for all the presented result 
we chose T = 0.3, which determines the shape of the PG 
density of states N P c(E) from Eq(4). 
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FIG. 1. (a) The normalized tunneling conductance 
i?jvdJ /dV as a function of the bias voltage in units of Ap with 
A s = 0.5 A P for different Z (= 0, 0.5, 1). N PG (E) (PDOS) 
is also plotted for reference; (b) The same as (a) but an- 
gle-averaged; (c) The same as (a) but with As = Ap; (d) 
The same as (c) but angle- averaged. 

In Fig. 1(a) we show the normalized tunneling conduc- 
tance R^dl /dV at zero temperature with the supercon- 
ducting gap As = 0.5Ap for different barrier potentials 
(Z=0,0.5, and 1). We also plot N PG (E) with T = 0.3 for 
comparison. The results can be trivially understood. As 
well known from the BTK theory, for small value of Z, we 
see the enhanced conductance below As due to the An- 
dreev reflection and this effect quickly disappears with in- 
creasing Z value. The main difference of our results from 
the conventional BTK calculation is that we modulate 
the conventional BTK tunneling conductance by multi- 
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plying with Npg{E). As a result the enhanced conduc- 
tance below Ag and the sharp peak structure at E = As 
appear inside the PG density of states (As < Ap). In 
Fig. 1(b) we show the same calculations but angle aver- 
aged to simulated the surface roughness in real tunneling 
experiment. The main features are the same as Fig. 1(a) 
and the line shapes become more rounded off. In partic- 
ular the result of Z=l case in Fig. 1(b) looks quite similar 
with the recent tunneling experiment by Krasnov et al. 

I 

In Fig. 1(c) and Fig. 1(d) we show the results of same 
calculations as Fig. 1(a) and Fig. 1(b), respectively but 
with As = Ap. When Ag has the same value as Ap 
the line shape of the tunneling conductance looks simi- 
lar to the conventional NIS junction, in particular with 
a large Z value. The effect of the PG is just to enhance 
the over-all line shape of the conductance. For smaller Z 
value the conductance is still enhanced below Ag due to 
the Andreev reflection. However even the enhanced con- 
ductance at E — > limit is far below than 2 (for a conven- 
tional NIS junction it approaches 2 for E < As as Z — > 
0) because of the reduced DOS of PG origin. Again the 
angle averaged results in Fig. 1(d) shows a more rounded 
off line shapes. The angle averaged PG density of states 
(PDOS) is shown for comparison in Fig. 1(b) and (d). 

In Fig2.(a-d), we plot the normalized conductance with 
varying A s (=0.5,l, and 1.5 A P ). When Z=l (Fig.2(a) 
and Fig. 2(b)) the results are easily understood. For 
Ag < Ap the distinct peak structure due to SC appears 
inside the PG as explained above. And when As > Ap 
there is no more distinct peak structure and the SC 
feature overwhelms the PG structure. When Z=0.1 
(Fig. 2(c) and Fig. 2(d)) the conductance line shapes look 
more peculiar, but basically it can be understood as an 
overlap of the Andreev enhanced conductance below As 
and the PG density of states Npo(E). G. Deutscher et 
al. H reported that there are two energy scales observed 
in tunneling experiments, and depending on the barrier 
strength (Z) only one of the two energy scales dominates 
the conductance line shape. In Fig2.(c) the line with 
As = 0.5Ap (solid line) shows a similar feature to the 
experimental data in Fig. 2 of Ref. m. 

In view of the experiments of Krasnov et al. M and 
G. Deutscher et al. ||, if our calculations have any rel- 
evance with underdoped HTSC, it would be the case of 
A 5 < Ap (Fig.l(a)(b) and Fig.2(c)(d)) and it will be 
interesting if more tunneling experiments with different 
Z parameters become available in near future. 



III. PG MODEL II:SDW 

In this section wc consider another possibility of the 
PG, namely, in which the correlation, which induces the 
PG, coexists with a SC correlation below T c and the PG 
correlation interplay with the SC correlation. In this case 



the key difference from the case (I) is that either As or 
Ap does not show up as separate features in the tunnel- 
ing conductance but the total gap A tota i 



A 2 
p 



shows up as a result of the interplay of two gaps. However 
more detailed line shape of the tunneling conductance re- 
veals the existence of two gaps, As and Ap. 




H ■ '. 
1 .' \ 




' 




I 1/ 

1 








//• 


(B) Z=1 


//-' 


(angle-ave) 







E/A_ 



E/A„ 



> 
■a 



■a 
c 

CC 2 





•' / y 




a v 




■ // 






(D) Z=0.1 




(angle-ave) 






1 


2 


E/A 


p 



E/A„ 



FIG. 2. (a) The normalized tunneling conductance 
RNdl/dV as a function of the bias voltage in units of Ap 
with Z=lfor different A s (= 0.5, 1, 1.5 A P ). iV PG (S)(PDOS) 
is also plotted for reference; (b) The same as (a) but an- 
gle-averaged; (c) The same as (a) but with Z=0.1; (d) The 
same as (c) but angle- averaged. 

Specifically, we assume the SDW correlation as the ori- 
gin of PG. The reason for this assumption is of course 
that there is a strong antiferromagnetic (AFM) correla- 
tion in underdoped HTSC compounds. Particularly for 
the underdoped HTSC compounds there is clear experi- 
mental evidence from neutron scattering fl5[ ] for the co- 
existence of the AFM correlation and superconductivity 
below T c and also the AFM correlation length becomes 
much larger than the SC correlation length. This enable 
us to treat the short range AFM correlation by SDW 
from the point of view of the SC correlation. 
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A. tunneling with SDW only 

As a prelude we first generalize the BTK theory to the 
NISDW (Normal metal-Insulator-SDW) junction. For 
simplicity we consider only a one dimension model and 
assume a commensurate SDW state with Q = kp = tt ja 
(a is the lattice distance) and also neglect the Fermi sur- 
face mismatch between the normal metal and the SDW 
state. 

The wave functions of the lefthand (normal metal) and 
the righthand (SDW) sides of the tunneling barrier is 
written as, 

qr L (x) = e ikx + Re~ ikx 

- T[a{k)e~ lkx + b(k)e~ l(k - Q)x ] (5) 

where a(k) 2 — ^jf and b(k) 2 = are the Bogoliubov 
coefficients of the SDW state, e = ^/£ 2 + A SDW , and 
By matching the boundary conditions 
of the wave functions at x = 0, we obtained R and T as 
follows. 



R(E) 
T{E) 



For e 2 > A| w , 



(a + b)iZ — a 
(a + b)iZ + b 
1 

(a + b)iZ + b 



(6) 
(7) 



e 2 < 



A SDW usin S a 



R\ 2 is easily calculated, and for 

> ^VgW^g and b 2 = 



2 c 



e+iy/A 2 SDW - 
2e 



-, we can show |i?| 2 = 1. 



The tunneling conductance of NISDW junction with 
different Z (=0,0.3, and 1) are shown in Fig.3(a,b,c) as 
dash-dot lines (As = case). The tunneling conduc- 
tance line shape is qualitatively different from the actual 
density of states Nsdw(E), which would have been ob- 
tained by the tunneling Hamiltonian method. Here we 
clearly observe that the BTK theory and the tunneling 
Hamiltonian method give different results for tunneling 
conductance. And it has been already known that these 
two methods describe different physical processes for the 
tunneling phenomena EJ. In particular, the tunneling 
Hamiltonian method, although physically appealing, has 
never been justified in regard with the tunneling trans- 
fer matrix element T (it is always taken as a constant). 
Then another question may arise: why then for the BCS 
state does the BTK theory give a qualitatively similar re- 
sult as the tunneling Hamiltonian methods at least for a 
large Z limit where the Andreev reflection process is sup- 
pressed, but not for the SDW state ? The reason is that 
for the BCS state the quasiparticle is a superposition of 
a momentum k particle and a momentum — k hole, and 
all together it carries the flux of the same momentum k, 
while for the SDW state the quasiparticle is a superposi- 
tion of momenta k and k + Q particles and therefore the 
quasiparticle in SDW state doesn't carry a single momen- 
tum. For the BTK theory of tunneling the main physics 



is the flux conservation described by Liouville's theorem, 
therefore the BTK theory sensitively traces the correla- 
tion of different momentum states in the SC or SDW 
states, while the tunneling Hamiltonian method is com- 
pletely blind of the momentum correlation. In the case 
of SC (BCS) state, we are in a fortune situation as de- 
scribed above, therefore momentum correlation doesn't 
play an important role but particle-hole correlation plays 
an important role in the BTK theory, which provides the 
main difference between the BTK theory and the tunnel- 
ing Hamiltonian method through the Andreev scattering. 
Now for the case of SDW state, the momentum correla- 
tion is important but not the particle-hole correlation; 
therefore it is expected that the two methods will give 
different results. 

An important question is then which theory should be 
trusted. As described above, the tunneling Hamiltonian 
method has never been justified, while there is no logical 
flow in the BTK theory based on the flux conservation. 
Therefore we think unless the tunneling interface is very 
rough (of course we need to estimate how much rough is 
rough) the BTK theory becomes more trustable. 



B. tunneling with SDW+SC 

Now we will consider the tunneling junction of a 
normal metal-insulator-SDW+SC (NISDWS). Again for 
simplicity we consider only a one dimensional model. In 
reality, as in HTSC, the analysis of tunneling in the two 
dimensional SDW state is more complicated, in particu- 
lar when SDW+SC is considered. However our one di- 
mensional model is enough to provide a qualitative un- 
derstanding of the interplay of SDW and SC correlations 
in tunneling process. 

Once the SDW state is formed, there are two branches 
of quasiparticles created: a+ t k = a(k)ck + b(k)ck+Q and 
ot— t k = b(k)ck — a(k)ck+Q, where a(k) and b(k) are the 
Bogoliubov coefficients as defined above. Now we as- 
sume that the superconducting pairing occurs between 
a+ ; fc and a+.-k, and also between and but 
not between a+^k and ct- t -k, for example. Although this 
assumption is mainly for simplicity of the analysis, it 
is known that the pairing interaction between different 
branches are much weaker in Hubbard model p7y . Also 
strictly speaking the SDW state with a commensurate 
wave vector (Q = ir/a) is an insulator with a fully devel- 
oped gap below Asdw = Ap, in contrast to the PG state 
in HTSC where residual density of states still remain and 
the system remains a metal. Therefore it is clear that 
our SDW state only mimic a state with a short range 
AFM correlation of the underdoped HTSC and it is jus- 
tified by the fact that the AFM correlation is effectively 
long range from the viewpoint of the SC correlation. A 
main drawback of this assumption is that the effect of the 
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residual density of states below PG in the tunneling con- 
ductance is completely missing. In summary, the purpose 
of this section is to study a tunneling characteristics of a 
SDW+SC state due to the interplay of two correlations, 
but not all the details of real materials. 

Now for the BTK theory the wave functions of the left- 
hand (normal metal) and righthand (SDW+SC) sides of 
the tunneling interface are written as, 



\A\ 2 = 



\R\ Z = 



A 2 
H ' 



(A P + 2Z y/A 2 - E 2 ) 2 + AZ 2 {Zy/A 2 -E 2 + A P ) 2 



H 



H = E 2 + [{2Z 2 + l)^A 2 -E 2 + 2ZA P } 2 . 



(12) 



■ikx 



ikx i a { rp\ ( \ ikx 



9 R {E)=C(E)fya + {-k) + D(E)(l) 



a_(-fe). (8) 



VZ 2 + &SDW + A S and £ 



where E 2 = ^fe 2 
h 2 k 2 /2m — fj,. All k w kp approximation is taken; the 
error of the approximation is O ( A to tai /Ep ) . Also a short- 
hand notation of the wave functions are 



a+(-k) = a{e)e~ lkx 
a-{-k) = b{e)e- tkx - 



a{ e )e-^ k - Q >. 



(9) 



As usual the SC Bogoliubov coefheients are u(E) 2 = 
(E + e)/2E and u(E) 2 = (E - e)/2E. By matching 
boundary conditions, which is a lengthy but straight for- 
ward calculation, we obtained the coefficients of A(E) 
and R(E) as follows. 

uv 

A = — , 
F ' 



R = 



(u 2 - v 2 )[Z 2 (b 2 - a 2 ) + iZ{b 2 - a 2 + 2ab) - ab] 



F=(u 2 - v 2 )[Z 2 (b 2 - a 2 ) + iZ ■ 2ab\ + b 2 u 2 + a 2 v 2 . (10) 



In order to calculate \A\ 2 and |i?| 2 we need to divide 
the region of E into three regions. For region (I), where 
E > A total (= ^A 2 p + A 2 ), e = ^E 2 - A 2 and ^ = 
For region (II), where A s < E < A totah 



^JE 2 - A 2 



total ' 



e = ^E 2 - A 2 and C = 



total 



E 2 , and finally 

for region (III), where E < A s , e = i^/A 2 - E 2 and 
But in the final results of \A\ and \R\, 



W^total-E 2 - 



e factors are all cancelled out, so that the region (II) and 
(III) are not distinguished. Therefore the tunneling char- 
acteristics shows distinct change only across E = A tota i 
but no change across either E — As or E = Ap. 

The final results of the Andreev (A) and normal reflec- 
tion (R) coefficients are: for region (I), 

A 2 



\R\ 2 = 



(2Z 2 y / E 2 - A 2 - A P ) 2 + AZ 2 {y/E 2 - A 2 + A P ) 2 



G 



G = \{2Z 2 + l)^E 2 - A 2 + Ef + 4Z 2 Ap; 
and for region (II) and (III), 



(11) 



> 

3 1 

■a 

c 

DC 









£ 


(A) Z=0 


i'. 


A a =2 




A s =1 


1 
i 


A B =0.5 


i 


— - A,=0 


) 1 


2 3 ' 



-^i 




..A 




i: 
i 


(B) Z=0.3 | 



E/A n 



2 

E/A„ 




FIG. 3. (a) The normalized tunneling conductance 
R N dI/dV for NISDWS (Normal metal-Insulator-SDW+SC) 
junction as a function of the bias voltage in units of Ap with 
Z=0 for different A s (=0, 0.5, 1, 2 Ap); (b) The same as (a) 
but with Z=0.3; (c) The same as (a) but with Z=l. 



In Fig.3(a-c) we plot the numerical calculations of the 
normalized tunneling conductance as a function of bias 
voltage (in unit of Ap) with varying size of As (=0,0.5,1, 
and 2) for different barrier potentials (Z=0, 0.3, and 1). 
The main features of the tunneling conductance are very 
different from a pure SC case. When As is much bigger 
than Ap (say, As/Ap — 2 in Fig. 3) the line shape is 
similar to the pure SC case. But even in this case the 
position of a gap in the conductance is determined by At 
and the presence of SC gap Ag only show up through an 
Andreev scattering coefficient, which enhances the con- 
ductance below the total gap. On the other hand when 
As < Ap, the line shape looks closer to the pure SDW 
case, which is qualitatively different from a real density 
of states as explained in the previous section. Again the 
presence of SC gap As show up only through an Andreev 
scattering and enhances the conductance below the total 
gap. 
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FIG. 4. (a) The normalized tunneling conductance 
R N dI/dV for NISDWS (Normal metal-Insulator-SDW+SC) 
junction as a function of the bias voltage in units of Ap with 
A s = 0.5A P for different Z (=0, 0.3, 1, 2); (b) The same as 
(a) but with As = Ap . 

In Fig. 4(a) (b) we plot the same calculations in different 
organization. Here we vary Z values with fixed As to see 
more clearly the effect of varying strength of the tunnel- 
ing barrier. Main features are already explained above. 
Both in Fig. 4(a) and (b), the normalized conductance 
below Atotal is enhanced by the Andreev reflection and 
it increases as Z — ► 0. When Z—0, we can have a simple 
results of \A\ 2 = §| and |i?| 2 = for E < A t 



A 2- 



*total ■ 



These expressions of \A\ 2 and |i?| 2 clearly show the con- 
trasting role of the SDW gap (normal reflection) and the 
SC gap (Andreev reflection); the normalized conductance 

is given by R N dI/dV{E < A t ) = ==£ in this limit. In- 
creasing Z from zero, it continues to suppress the Andreev 
reflection and enhances the normal reflection. 

In view of current tunneling experiments, if our results 
of the NISDWS junction has any relevance with the un- 
derdoped HTSC, only possibility is the case of As > Ap 
at low temperature. In order to realize this possibility 
starting with a sizable Ap and a zero As above T c , we 
have to imagine that Ap gradually decreases and As 
gradually increases as temperature is lowered. Because 
this scenario is quite unlikely from the currently available 
experiments, we can conclude either that the PG is not 
caused by an AFM correlation, or that even if the origin 
of PG is an AFM correlation the AFM correlation is not 
strong enough to have any significant interplay with a SC 
correlation. 



IV. CONCLUSION 

We considered two phenomenological models of the PG 
state in underdoped HTSC and studied the characteris- 
tics of tunneling conductance at zero temperature by gen- 
eralizing the BTK theory. In model I, we assumed that 



the PG is a simple suppression of density of states of a 
unknown origin and when the system goes to SC state 
there is no direct interplay between the PG correlation 
and the SC correlation. In this case the characteristics 
of tunneling conductance is a simple superposition of a 
standard tunneling conductance with SC gap As and the 
PG density of states Npq(E). Despite the simplicity of 
the model, the results seem to explain some of the re- 
cent tunneling experiments by Krasnov et al. Q and G. 
Deutscher S, which indicate non-superconducting ori- 
gin of the PG. In the model II, we assumed that the PG 
is caused by an AFM correlation and it is simulated by 
SDW state. Below T c the SDW correlation and SC cor- 
relation show an interesting interplay. As a result, the 
tunneling gap is given by A tota i — y/Ag + Ap and in- 
dividual gaps, As and Ap, do not show up explicitly. 
In particular when Ap^sdw > As, the line shape of the 
tunneling conductance looks qualitatively different from 
a conventional NIS junction. In view of available tunnel- 
ing experiments in underdoped HTSC, the relevance of 
the model II is possible only when Apsrw < As below 
T c , which is quite unlikely at present. To clarify the issue 
of the PG more experiments on tunneling are essential 
and our study should serve as a useful benchmark. 
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